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ABSTRACT 

The polar disk galaxy NGC 4650A was observed during the commissioning of the Multi Unit Spectroscopic Explorer (MUSE) at 
the ESO Very Large Telescope to obtain the first 2D map of the velocity and velocity dispersion for both stars and gas. The new 
MUSE data allow the analysis of the structure and kinematics towards the central regions of NGC 4650A, where the two components 
co-exist. These regions were unexplored by the previous long-slit literature data available for this galaxy. The stellar velocity field 
shows that there are two main directions of rotation, one along the host galaxy major axis (PA. = 67 degrees) and the other along the 
polar disk (PA. = 160 degrees). The host galaxy has, on average, the typical pattern of a rotating disk, with receding velocities on the 
SW side and approaching velocities on the NE side, and a velocity dispersion that remains constant at all radii (cTstar ~ 50 - 60 km/s). 
The polar disk shows a large amount of differential rotation from the centre up to the outer regions, reaching V ~ 100- 120 km/s 
at R ~ 75 arcsec ~ 16 kpc. Inside the host galaxy, a velocity gradient is measured along the photometric minor axis. Close to the 
centre, for R < 2 arcsec the velocity profile of the gas suggests a decoupled component and the velocity dispersion increases up to 
~ 110 km/s, while at larger distances it remains almost constant (o-gas ~ 30 - 40 km/s). The extended view of NGC 4650A given by 
the MUSE data is a galaxy made of two perpendicular disks that remain distinct and drive the kinematics right into the very centre 
of this object. In order to match this observed structure for NGC 4650A, we constructed a multicomponent mass model made by the 
combined projection of two disks. By comparing the observations with the 2D kinematics derived from the model, we found that the 
modelled mass distribution in these two disks can, on average, account for the complex kinematics revealed by the MUSE data, also in 
the central regions of the galaxy where the two components coexist. This result is a strong constraint on the dynamics and formation 
history of this galaxy; it further supports the idea that polar disk galaxies like NGC 4650A were formed through the accretion of 
material that has different angular momentum. 

Key words. Galaxies: kinematics and dynamics; Galaxies: individual: NGC 4650A 


1. Introduction 


The polar disk galaxy NGC 4650A is the prototype of the class of 
polar-ring galaxies (PRG, see |Iodice|2014[ for a review). Since 
it is characterized by two components, the central host galaxy 
(HG) and the polar structure that rotate on almost perpendicular 
planes, it is also classified as a multispin galaxy ( |Rubin|1994] ), a 
class that includes all systems with a kinematically distinct com¬ 
ponent of gas and/or stars with a variety of inclination angles 
and radial extent with respect to the host galaxy. Two main pro¬ 
cesses were proposed for the formation of the decoupled com¬ 
ponent in multispin galaxies: accretion of material (gas and/or 
stars) from outside (from other galaxies and/or environments) 
by the pre-existing host galaxy, or merging galaxy (see |Comb^ 
|2014|[Iodice|2014[ for reviews). In the cold dark matter scenario 
for galaxy formation, such gravitational interactions play a fun- 


* This work is based on observations taken at the ESO La Silla 
Paranal Observatory within the MUSE Commissioning. 


damental role in defining the morphology of “normal” galax¬ 
ies, in particular in the building up of spheroids. Thus, in this 
framework, the study of multispin galaxies, both at low and high 
redshift, can shed light on the main processes at work during 
galaxy interactions and on the influence of the environment (see 
|Conselice|2014[ for a review). Moreover, the existence of two or¬ 
thogonal components of the angular momentum makes the PRGs 
the ideal laboratory in which to derive the 3D shape of the dark 
matter spheroid that dominates the gravitational potential (see 
|Arnaboldi|2014[ for a review). 

In the last ten years, a large amount of data (images in the op¬ 
tical and near-infrared bands, emission and absorption line long- 
slit spectroscopy and HI radio emission) has been collected for 
NGC 4650A. In the combined analysis of the whole data set, 
by studying morphology, kinematics, and stellar population, the 
structure and formation for this fascinating galaxy was outlined 
very accurately. In detail, NGC 4650A is the first PRG classi¬ 
fied as a polar disk galaxy; in these galaxies stars and dust in 
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the polar structure are traced inward within the HG, down to 2. Observations and data reduction 


the galaxy centre (Arnaboldi et al. 1997 

lodice et al. [2002a ; 

[Gallagher et al. 2002[ Swaters & Rubin 

2003[). Thus, the two Observations of NGC 4650A were taken with the Multi Unit 

components co-exist in the inner regions." 

rhe main properties of Spectroscopic Explorer (MUSE) mounted on Unit Telescope 4 


the polar structure are i) the large amount of HI gas, 10^^ M©, 
which is four times more extended than the optical counterpart 
(i.e. up to 40 kpc from the centre); ii) the differential rotation; 
and in) the subsolar metallicity Z = 0.2Zo, which is constant 
over the whole extension of the disk (Spavone et al. 2010| ). Given 


its spheroidal shape, the HG in NGC 4650A has been classified 
as a SO-like system. This morphological classification was dis¬ 
puted by |Iodice et al.| ( |200^ on the basis of the stellar kinematics 
along the spheroid’s main axes. In fact, even if this component is 
rotationally supported, with a maximum rotation velocity along 
the major axis of V ^ 80 - 100 km/s, the velocity dispersion 
remains almost constant (cr ~ 65 km/s) at all radii and along 
both axes. These measurements put the HG far from the Faber- 
Jackson rela tion for early-type galaxies (see |Iodice et ar]|2006[ 
|Iodice||2014| ), which shows that NGC 4650A has lower central 
velocity dispersion than that measured for spheroids of compa¬ 
rable luminosity. 

The structure, kinematics, and metallicity derived for 
NGC 4650A are consistent with the formation of a polar disk 
galaxy through the accretion of gas from outside, either from a 
cosmic filament or from a gas-rich donor with a different angular 
momentum ( [Spavone et al.|2010t|Combes|2014| ). 


(UT4) of the ESO VLT at Paranal, Chile, during the instrument 
commissioning run on 8 and 12 February 2014. Observations 
were acquired in the wide field mode, having a 1 x 1 arcmin 
field of view and a scale of 0.2 arcsec/pixel. The spectral range 
is 4800-9 300A and the instr umental dispersion is about 55 km/s 
at 5200A ( Bacon et al.|2010 ). 


From dynamical studies for the kinematics of NGC 4650A, 
the best models predict a flattened E6-E7 dark halo, wit h its ma- 
jor axis aligned along the plane of the polar disk itself ( [Combes j 
' & Arnaboldi|1996[[Napolitano et al.|2014[ ). The biggest uncer¬ 


tainties in these mass models reside i) in the limited spatial cov¬ 
erage of the long-slit spectra and ii) in the contamination of the 
HG kinematics by the stars in the polar disk, in the regions where 
the two components coexist. The long-slit spectra cannot trace 
both the S-shape morphology of the polar disk towards the cen¬ 
tre of the g alaxy ([lodice et ar][20Q2b[ ) and the warped arms at 
large radii ( [S waters & Rubin[[2003[ ), thus the star and gas mo¬ 
tions for this component are not mapped in these regions. 

Recently, [Coccato et al.[ ( [201^ made a first attempt to disen¬ 
tangle the contribution of the stars in the two components (HG 
and polar disk) along the HG minor axis. The main results indi¬ 
cate that along this direction i) the significant rotation measured 
seems to be an intrinsic property of the HG rather then an artefact 
caused by the contamination of the stars in the polar disk and ii) 
the polar disk shows a kinematic decoupling in the inner regions, 
which is counter-rotating with respect to the outer regions and 
the HG. These findings pointed out again the complex structure 
of NGC 4650A in the regions where the two components coexist 
and suggested that the determination of the internal structure of 
NGC 4650A clearly requires integral-field spectroscopy. 

In this work we present new observations taken for 
NGC 4650A with the Multi Unit Spectroscopic Explorer 
(MUSE) at the ESO Very Large Telescope (VLT). The wide 
wavelength range and the high spectral resolution combined with 


the large field of view and fine angular sampling of MUSE ( [Ba- 
[con et al.|2010[ ) allow us to derive the first 2D field of the rota¬ 
tion velocity and velocity distribution for both stars and gas in 
NGC 4650A. Thus, we are able to analyse the entire 2D struc¬ 
ture (morphology and kinematics) of both HG and polar disk, in 
particular in the regions close to the galaxy centre where both 
components contribute to the light distribution and kinematics. 

In this paper, we adopt a distance to NGC 4650A of 
43.4 Mpc and Hq = 73 km s“l Mpc“L This implies that 
1 arcsec= 0.21 kpc. 


The data consist of 12 dithered exposures of 10 minutes each, 
divided into three pointings. The central pointing, obtained by 
combining six exposures, covers the host galaxy of NGC 4650A 
and part of the polar structure. The northern and southern arms of 
the polar disk are covered by combining three exposures for each 
side. A standard calibration was adopted for this run, which in¬ 
cludes bias, fiat field, and comparison arcs. A standard star (CD 
71) was observed at the beginning of each night to correct for 
instrument response across the wavelength range. The average 
seeing was FWHM ~ 0.6 arcsec. 

Data reduction was performed using the MUSE pipeline ver¬ 
sion 0 .18.5 ([Weilbacher et al.[2014T ) under the REFLEX environ¬ 
ment ( [Freudling et al.[2013 ). The sky background was evaluated 
on the regions of the field of view where the galaxy contribu¬ 
tion was negligible. The different pointings were aligned using 
foreground stars as reference, and then combined to create a fi¬ 
nal datacube. The alignment consisted of offsets along the RA 
and DEC direction, and a ~ 3 degree rotatioi£] We noticed that 
imperfect sky subtraction left strong residuals for wavelengths 
larger than -7000 A. We therefore concentrated our analysis on 
the wavelength range 4780 A - 6850 A. 

The final MUSE mosaic of NGC 4650A is 1.5 x 2.5 arcmin 
(see Fig. [T]). The MUSE optical r-band image of NGC 4650A 
(see left panel Fig. shows clearly that the light distribution is 
dominated by two almost perpendicular disk-like components. 

The MUSE map of the Ha emission is shown in the right 
panel in Fig. [T] As already known for this galaxy ( [Arnaboldi| 
[et al.|1997} Swaters & Rubin|2003] ), the neutral and ionized gas 
component resides completely in the polar disk. The MUSE data 
reveal the whole 2D distribution of the ionized gas and that the 
emission from the star forming regions is also quite strong close 
to the centre of the galaxy. 

The MUSE data for NGC 4650A are consistent with the view 
given by previous HI data ( [Arnaboldi et al.[[1997[) a nd the op- 
tical and n ear-infrared images ( [lodice et ^ 2002a[ [Gallagher 
[et al.|2002| ). They showed that the S-shape observed in the near- 
infrared images, the dust and bright blue knots detected in the 
optical data, and the HI iso-velocity contours are reconciled by 
the presence of spiral arms in the almost edge-on polar disk 
( [Arnaboldi et al.[[r997[ ). These features originate near the cen¬ 
tre of the galaxy, the southern arm passes in front of the HG and 
its light and dust affect the underlying HG component on the SW 
side. The structure of the spiral arms derived by the best fit of the 
HI data is plotted on the MUSE r-band image of NGC 4650A 
(see left panel in Fig.[T]). The Ha emission shows several bright 
knots in correspondence with the outer parts of the NW and SE 
arms. 


^ Because the pipeline currently does not handle rotation when com¬ 
bining multiple frames, it was necessary to overwrite the header key¬ 
words on the raw frames accordingly. 
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lodice E. et al.: Two perpendicular disks in NGC 4650A 
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Fig. 1. Left panel - MUSE mosaic of NGC 4650A in the SDSS r-band. The mosaic covers an area of 1.5' x 2.5'. The X-axis is the Right Ascension 
(hr min sec) and the Y-axis is the Declination (in degrees). Dashed boxes indicate the single MUSE pointings. Dashed and dotted arcs trace the 
structure of the spiral arms (see text). Right panel - Her emission from MUSE data. North is up and east is on the left. 


3. Two-dimensional kinematics of stars and gas 

3.1. Measurements 


To measure the absorption line kinematics, it was necessary to 
add the spectra of adjacent spatial pixels in the field of view. To 
this aim, we adopted the Voronoi Binning Scheme, using the IDL 
implementation by |Cappellari & Copin| ( |2003] ). The ionized gas 
emission lines were strong enough to be measured on individual 
spatial pixels, except for the innermost 5", where we used the 
spatial bins defined in the Voronoi Tesselation. 

We first measured the stellar and ionized gas kinematics on 
the sp atial bins defined by the Voron oi Tesselatio n, using the 
pPXF ( [Cappellari & Emsellem||2004| ) and Gandalf ( |Sarzi et aL 


|2006| ) IDL fitting routines. We then used the measured stellar 
kinematics to constrain the stellar continuum in each single spa¬ 
tial pixel and subsequently measured the emission line kinemat¬ 
ics. 


In both the absorption and emission line fits, we include mul¬ 
tiplicative polynomials of degree 8 to account for the shape of 
the stellar continuum. In the pPXF fitting, we used the library of 
single stellar population spectra from |Vazdekis et al^ ( |20I2| ). 

Each emission line was fitted independently; the mean 
ionized gas kinematics was derived by averaging the re¬ 
sults for each emission line. The fitted emission lines are 


HyST4861, [OIII]TT4959,5007, [NI]TT5198,5200, [HeI]T5876, 
[NaI]TT5890,5896, [OII]TT6300,6364, [NII]T6548, HaT6563, 
[NII]T6583, [SII]T6717, and [SII]T6731. The amplitude versus 
noise of each line was used as a weight in the average; we did 
not consider emission lines whose amplitude-to-noise ratio was 
lower than 4. The mean velocity and velocity dispersion are the 
luminosity weighted mean of the lines that satisfy the above cri¬ 
teria. For the majority of the spatial bins, only the Ha, Hj3, [Nil], 
and [OIII] lines satisfy the above criteria. 


3.2. Stellar Kinematics 


The 2D maps of the line-of-sight velocity and velocity disper¬ 
sion for stars in NGC 4650A are shown in Fig.|^ The two main 
directions of the star rotation are along the HG major axis (where 
the position angle is P.A. = 67 degrees) and along the polar disk 
(P.A. = 160 degrees), which is almost perpendicular to the HG 
(see left panel in Fig.[^. 

Inside the HG, the velocity dispersion remains quite con¬ 
stant at all distances from the galaxy centre and P.A.s, at cr ~ 
50 - 60 km/s, while it is larger along the polar disk up to 
(T ~ 100 km/s in the outer regions (see right panel in Fig.[^. The 
increase in the stellar velocity dispersion along the polar disk di¬ 
rection was already found from the long-slit data by [lodice et al. 
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( [2006| ). By using this data set, after separating the contribution 
by the stars in the HG from that of the stars in the polar disk, 
[Coccato et al.| ( |2014| ) concluded that the increase of velocity dis¬ 
persion observed along the polar disk is an intrinsic property. 

Below, we focus on the description of the 2D kinematics in 
the regions where the two components, HG and polar disk, co¬ 
exist. In particular, we discuss how features and perturbations in 
the velocity field of stars can be reconciled by the structure of 
the light distribution at the same radii. 

HG major axis - The overall iso-velocity contours inside the 
HG resemble the typical pattern of a rotating disk, with receding 
velocities on the SW side and approaching velocities on the NE 
side. Nevertheless, the MUSE measurements show that along the 
HG major axis the iso-velocity contours on the side indicates a 
flatter rotation curve than the corresponding regions on the east 
side (see top left panel in Eig.|^. Along this direction, the light 
distribution is also not symmetric with respect to the centre: the 
optical image shows that at 7? > 10 arcsec the SW isophotes are 
rounder than those at the same distance from centre on the other 
side, which appear more disky (see left panel in Eig. [^. This 
asymmetry in the light is more evident from the high-frequency 
residual imag^of the MUSE r-band mosaic where the disk of 
the HG is more extended on the NE side (R ~ 20 arcsec) than 
on the SW side (R < 16 arcsec, see the top left panel in Eig.[^. 
This is a strong indication that the asymmetry in the kinematics 
reflects an asymmetry in the structure of the HG disk. This effect 
can be quantified by looking at the folded light profiles and the 
folded line-of-sight velocity profiles extracted along the HG ma¬ 
jor axis, shown in Eig.|^ As noticed in the velocity field, the NE 
velocity profile rises at all radii, while the SW velocities tend to 
have constant values forR > 10 arcsec. The light profiles have 
a different slope forR > 16 arcsec, where the NE profile has a 
steeper decline than that the SW profile. 

Einally, we found that the kinematic and photometric major 
axis are not coincident. We measured the kinematic P.A. wofile 
by using the software Kinemetrj^ and it is shown in Eig.[^ The 
average value in the range of radii inside the regions of the HG 
that are not perturbed by the polar disk (4 < /? < 12 arcsec) is 
P.A.kin = 72.2 + 0.6 degrees. This was compared with the pho¬ 
tometric P. A. profile derived by the fit of the isophote performed 
on the MUSE mosaic of NGC 4650A in the r-band. In the same 
range of radii the photometric PA. is about 5 degrees smaller 
than the kinematic P.A., being P.A.pho = 67 + 2 degrees (see 
Fig.0. 

HG minor axis - Another interesting feature in the stellar ve¬ 
locity field inside the HG is the non-zero velocity curve along 
the PA.=157 degrees, which corresponds to the HG minor axis 
and to the polar disk major axis (see top left panel in Pig. [^. 
This suggests some rotation along this direction that was already 
found in the velocity curves derived from the long-slit data. By 
separating the contribution of the velocity of stars in the polar 
disk from that in the HG, which co-exist along this direction, 
Coccato et al.| ( |2014| ) suggested that this could be a real addi¬ 


tional feature of the galaxy structure. 

Perturbations due to the polar disk - The co-existence of the 
polar disk in the inner regions of the HG (as shown by the high- 
frequency residual image in the top right panel in Pig.[^ gener¬ 
ates some perturbations to the velocity distribution of stars. 


^ This is the ratio of the original reduced image with a smoothed one, 
where each original pixel value is replaced with the median value in a 
rectangular window of 25 x 25 pixels. We used the IRAF task FME- 
DIAN to smooth the original reduced image. 

^ Kinemetry is a generalization of pho tometry to the higher m oments 
of the line-of-sight velocity distribution ( [Krajnovic et al.|2006| . 
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Along the HG major axis, for 3 < R < 6 arcsec, the iso¬ 
velocity contours are affected by the stars in the polar disk arms 
which pass in front of the HG on the SW side and behind the HG 
on the NE side (see upper panel in Pig. [^. The same features 
are also evident in the folded light profiles and rotation curves, 
shown in Pig.|^ 

On the NW side, at 7? ~ 7 arcsec from the centre (see top left 
panel in Pig.[^, the velocity field shows a region where the stel¬ 
lar velocity is blue-shifted with respect to those of the surround¬ 
ing area, which are red-shifted. In the same region, the velocity 
dispersion is larger than the average value at the nearby points 
(see right panel in Pig.|^. We checked that these measurements 
are not the result of a bad fit. They correspond to the region of 
the NW arm of the polar disk that emerges from the galaxy cen¬ 
tre on the north side and then goes down, passing in front of the 
HG (see top right panel in Pig.|^. The measured velocity and ve¬ 
locity dispersion in this spatial bin, where the contribution from 
the stars in the HG is negligible, have similar values to those in 
the southern arm of the polar disk at about 10 arcsec from the 
centre. 

The same explanation can be valid for another “out-of- 
bound” region on the SE side of the HG at 7? ^ 15 arcsec from 
the centre, where velocity is lower and dispersion is higher with 
than in the nearby regions (see top left panel in Pig. and the 
right panel in Pig.[^. This area can match the kinematics of the 
fainter southern arm of the polar disk which goes up and inter¬ 
sects the outer regions of the HG on the east side (see left panel 
in Pig. and top right panel in Pig.[^. In this case the velocity 
and velocity dispersion are also consistent with the values mea¬ 
sured for the SE arm of the polar disk. 


3.3. Comparison with the kinematics from long-slit data 


Prom the 2D kinematics of the stars in NGC 4650A we derived 
the profiles of the line-of-sight velocity and velocity dispersion, 
along the major and minor axis of the HG. These are shown in 
Pig. and are compared with the same quantities for this com¬ 
ponent in NGC 4650A obtained by the long-slit data take n with 
PORS2 at VET ( [lodice et al.|2006| ) and data published by jSack- 
|ett et al.| ( |1994| ). The agreement among the new and literature 


line-of-sight velo city profiles is very g ood. The velocity disper¬ 
sion measured by [lodice et ah] ( |2006| ) from the long-slit data is 
systematically larger for measurements based on the the absorp¬ 
tion lines around 5 lOOA. We understand that this is a systematic 
effect caused by the contamination of the Paschen lines in the 
calcium triplet (CaT) wavelength rang e, which leads to a h igher 
estimation of the velocity dispersion podice et al.||2006| ). The 
MUSE 2D map of the Ha emission also shows a strong intensity 
in the regions of the HG (see Pig.[^ right panel). This suggests 
the presence of hot and young stars and hence the contamination 
of the Paschen lines as the origin to the systematic effect that 
caused lager velocity dispersion measurements for the CaT lines 
versus those measured with the Mgb and Pe lines. The velocity 
dispersion values measured from the MUSE spectra are derived 
from the Mgb and Pe lines and they are smaller. 


3.4. Gas Kinematics 

The 2D maps of the line-of-sight velocity and velocity disper¬ 
sion derived by the emission lines in NGC 4650A are shown in 
Pig. [ 7 ] and in the bottom right panel in Pig. The ionized gas 
kinematics is mapped out to 75 arcsec from the centre, i.e. up to 
about 16 kpc, while the previous rotation curves from long-slit 
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Fig. 2. Map of the line-of-sight rotation velocity (left panel) and velocity dispersion (right panel) of stars in NGC 4650A. The cross marks the 
centre of the galaxy where the systemic velocity is V = 2875 + 3 km/s and velocity dispersion is cr = 60 + 6 km/s. 


data extend up to 10 - 12 kpc ( |S waters & Rubin|2003| ). The 2D 
velocity map of the ionized gas show that the polar disk has a 
large amount of differential rotation from the centre out to the 
outer regions, reaching V ~ 100- 120 km/s at R ~ 75 arcsec 
~ 16 kpc. This is consistent with the HI distribution that was 
well fitted by a model of disk in differential rotation (|Amaboldi| 
|etal.|1997| ). 


From the MUSE data of the ionized gas, the velocity of the 
approaching and receding arms on the NW and SE sides of the 
polar disk is measured for the first time. The gas velocity disper¬ 
sion remains almost constant (cTgas ~ 30-40 km/s) at all radii 
(see right panel in Fig. [ 7 ]). 

Along the polar disk, there are some regions where the line- 
of-sight velocity and velocity dispersion have different values to 
those measured in the nearby regions. This is particular evident 
on the east side and on the outer north-west side, where veloci¬ 
ties are lower (close to the systemic value) than in the adjacent 
regions. These features reflect the structure of the polar disk. As 
discussed in Sect.|^ the polar disk contains two prominent spiral 
arms that originate near the centre of the galaxy, and two fainter 
arms that are visible at larger radii. The fainter southern arm goes 
up on the east side and it overlaps the north arm at about 30 arc- 
sec from the centre (see Fig. and also left panel in Fig. [T]). 
As a consequence, in these regions the velocity of the gas has 
different values from that in the surrounding area. Consistently, 
the velocity dispersion also has different values to those in the 
adjacent regions. 

We measured the rotation curve and velocity dispersion pro¬ 
file along the polar disk major axis, at PA. = 160 degrees, and 
they are shown in the left panels in Fig. The line-of-sight ve¬ 
locity increases with distance from the galaxy centre. It reaches 
~ 100 km/s at ~ 50 arcsec on the NW side and ~ 120 km/s at 
~ 58 arcsec on the SE side. The rotation curve along the polar 
disk derived from the MUSE data is consistent with that derived 
by [Swaters & Rubin] ( [2003| ) from the [OIII] and Ha emission 


lines by using slit spectra. The above values are consistent with 
the HI velocities in the range 4 < R < 16 kpc|^ As pointed out 
previously, the velocity dispersion along the polar disk is almost 
constant (cr ~ 30 - 40 km/s) at all radii, except for the centre. 

3.5. Gas kinematics in the galaxy centre 

The MUSE data for NGC 4650A allow us to map the 2D kine¬ 
matics of the ionized gas all the way to the centre of this galaxy. 
Previous long-slit data traced the ionized gas kinematics only 
along the two arms of the polar disk where the line emission is 
stronger and so did not cross the galaxy centre. 

The MUSE data show that inside 5 arcsec (see bottom right 
panel in Fig. [^, the iso-velocity contours are almost parallel to 
the minor axis of the polar disk (i.e. P.A.=70 degrees), suggest¬ 
ing a cylindrical rotation inside this region. Deviations from such 
behaviour are observed closer to the centre, forR< 1 arcsec, 
where the iso-velocity contours show a modest twist. As noted 
for the velocity distribution of the stars (see Sect. |3.2| ), the ob¬ 
served distortions in the velocity field of the ionized gas are due 
to the presence of the several star forming regions along the po¬ 
lar disk, which are particularly intense when they are close to 
the galaxy centre (see right panel in Fig. and top right panel in 
Fig.|3l>. 

Inside 2 arcsec from the galaxy centre, the velocity dis¬ 
persion profile shows a steep increase up to a value of about 
110 km/s (see upper right panel in Fig. [^. In the same range of 
radii, a hint of a decoupled component is also observed in the 
velocity curve (see bottom right panel in Fig.[^. However, tak¬ 
ing into account that the fit is based on a few pixels and that the 


^ The velocities derived b y the HI inside 4 kpc from the centre are 
affected by large errors (see |Iodice et al.|2008] ), thus values can be com¬ 
pared with velocities from the MUSE ionized gas only in the range 
4 < R < 16 kpc. 


Article number, page 5 of 


14 




























A&A proofs: manuscript no. n4650a_MUSE_astroph 



2780 2800 2820 2840 2860 2880 2900 2920 2940 2960 i 



•100 -50 0 50 100 2750 2000 2050 2900 2950 3000 


Fig. 3. Top left panel. Enlarged image of the 2D map of the line-of-sight velocity of stars in the region of the HG. The iso-velocity contours are 
overlaid on the image, the minimum and maximum values are 2785 km/s and 2960 km/s, with a step of about 6 km/s. The cross marks the centre 
of the galaxy where the systemic velocity is V = 2875 + 3 km/s. The dashed lines indicate the photometric major and minor axes of the HG. Top 
right panel High-frequency residual image obtained by the ratio of the original reduced r-band MUSE image with a smoothed one, where each 
original pixel value is replaced with the median value in a rectangular window of 25 x 25 pixels. Bottom left panel Enlarged image of the 2D map 
of the velocity field derived by the model. The dashed line indicates the direction of the HG minor axis P.A. - 157°. The iso-velocity contours 
are overlaid on the image, the minimum and maximum values are -145 km/s and 145 km/s, with a step of about 13 km/s. See Sect.|^for details. 
Bottom right panel Enlarged image of the 2D map of the line-of-sight velocity of gas in the regions close to centre of the galaxy. The iso-velocity 
contours are overlaid on the image, the minimum and maximum values are 2800 km/s and 2920 km/s, with a step of about 6 km/s. In all images 
the cross marks the centre of the galaxy and north is up, east on the left. 


error bars are large, we cannot be certain that this feature cor¬ 
responds to a real component. One possible explanation of the 
higher velocity dispersion is the presence in the centre of two 
components that lead to a broadening of the emission line. By 
looking at the spectra in this region we did not find any sign of a 
double peak. Thus, by summing two Gaussians convolved with 
the MUSE resolution, we checked that, if two components exist, 
they should have a separation in velocity of V < 60 km/s; other¬ 
wise, for larger values a double peak should be observed in the 
emission line. Alternatively, the existence of an outflow along 
the line of sight could generate an increasing velocity dispersion 
of the gas towards the nucleus, and an almost zero velocity, as 
observed for NGC 4650A. Finally, the same effect could also be 
due due to chaotic motions of the gas in the nuclear regions. 


4. The multicomponent mass model for NGC 4650A 

In this section we present the numerical model made to repro¬ 
duce the observed structure and the 2D kinematics of stars and 
gas showed by the MUSE data for NGC 4650A. 

Photometry and spectroscopy (including new MUSE data) 
have shown that the light distribution and kinematics are con¬ 
sistent with two almost perpendicular disks where the two main 
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directions of the star rotation are along the equatorial host disk 
{P.A. = 67 degrees) and along the polar disk {P.A. = 160 de¬ 
grees). The two components co-exist in the central regions of 
the galaxy since stars and dust in the polar disk can be traced up 
to the galaxy centre, following a spiral pattern (see top left panel 
in Fig. 1^. 

Even if the morphology of the central HG resembles that of 
an SO galaxy (see left panel in Fig.[T]), the observed light distri¬ 
bution podice et al.|2002b) and the constant velocity dispersion 
of stars at all radii (as revealed by FORS2 and the current MUSE 
absorption line spectra) strongly suggested that this component 
can be considered a puff-up disk, rather than a hot spheroid. The 
disk nature of the HG is also evident from the high-frequency 
structure along the equatorial direction found in the residual im¬ 
age shown in the top left panel in Fig.[^ 

Based on these observational facts and taking into account 
that the contribution to the total velocity comes from stars in 
the HG and in the polar structure, we built a minimum model 
made of two perpendicular disks. The low value measured for 
the velocity dispersion of both components further supports this 
idea. 

The model does not include any rotation along the minor axis 
of the HG, the spiral arms in the polar disk, and a decoupled 
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Fig. 4. Folded line-of-sight velocity profiles of stars (left panel) and surface brightness profiles (right panel) along the HG photometric major axis 
{PA. = 67 + 2 degrees) from the MUSE data. 
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Fig. 6. Stellar kinematics along the major axis (left panel) and minor axis (rig ht panel) of the HG i n NGC 4650A derived by MUSE data (black 
points) in the Mgb wavelengt h region compared with the profiles obtained by |Sackett et al.| ( p~994| ), blue points, in the same region, and longslit 
spectra by [lodice et al.| ( [2006] ), red points, in the CaT wavelength region. 


component in the centre, as well as the misalignment between 
kinematic and photometric major axis of the HG, also revealed 
by the MUSE data. This approach allows us to verify whether 
the minor axis rotation in the HG and/or the decoupled central 
component are due to the projection of the orbits in the polar 
disk onto the HG ones. If this is the case, we should also mea¬ 
sure these features in the velocity and velocity dispersion fields 
derived by the model. Alternatively, they are to be considered 
intrinsic properties of the HG in NGC 4650A. 


To derive the velocity field in both equatorial and polar 
planes, one technique is to build a 3D gravitational potential 
by combining density models of the host galaxy and the polar 
disk, and then to search for the symmetric close orbits in both 
the equatorial and perpendicular planes by shooting test parti- 
cles (e.g. jSackett & Sparke||1990t jSackett et al.||1994t [Combes 

[& Arnaboldi[[1996[ ). This procedure is quite long, since parti¬ 
cles are launched tangentially on the major axis at each radius 
of the minor axis, and the launching velocity is adjusted until 
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Fig. 5. Kinematic (blue points) and photometric (black points) RA. 
profiles. The straight lines (blue and black) indicate the average values 
of the two P.A.s in the range of radii inside the regions of the HG that 
are not perturbed by the polar disk (4 < R < 12 arcsec). 


the closed orbit is found. Instead, we adopt here the epicyclic 
approximation in each plane (equatorial and polar) to compute 
the closed orbits in a non-axisymmetric potential (e.g. |Binney &| 
Tremaine|1987| ). Indeed, in their respective planes, the existence 
of the perpendicular component is equivalent to a weak bar, with 
zero pattern speed. 

In this work, we choose the multicomponent mass model 
previously developed by [Combes & Arnaboldi| \\996) . This 
model was proposed in order to revisit the 3D shape of the dark 
matter halo in NGC 4650A, and they found that the luminous 
mass closely fitted the kinematics along the major axis of the 
HG galaxy and the optical rotation curve of the polar disk, while 
to account for the high polar velocities at large radii given by the 
HI data, a dark matter halo is required from 8-10 kpc. The main 
result of this study was on the shape of the dark matter halo: the 
best fit to data with the least amount of dark matter was obtained 
with a flattened E6-E7 dark halo, having its major axis aligned 
along the polar disk major axis. 

In this work we adopt the same configuration for 
NGC 4650A as in the model proposed by [Combes & Arnaboldi| 
( |1996[ ), i.e. two luminous components HG and a polar disk plus a 
flattened dark matter halo, with only slight modifications to take 
into account the extended and deeper photometry in the I and K 
bands published by [Iodice et al.[ ( [20'02b[ ). 

As in [Combes & Amaboldi[ ( [1996[ ), according to the light 
distribution, the central host galaxy can be described as a very 
small bulge plus a thick disk. The bulge is very light, and its 
influence might be ignored. It is represented here by a Plummer 
component, having a characteristic radius vb and a total mass 
Mb. The host disk is represented by a double exponential, with 
mass Mj, characteristic scale r^, and height h^. Ah the above 
parameters are given in Table 

The stellar and gaseous polar structure are represented by 
different Miyamoto-Nagai potential-density pairs ( [Miyamoto &[ 
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[Nagai|1975] ). Table indicates their masses Mpo, their common 
height h, and their two internal and external radii ri and r 2 . The 
model of the polar disk does not account for the observed spi¬ 
ral arms and for the several star forming regions that strongly 
perturb this component. 

A dark matter halo is added to account for the high velocities 
in the outer parts. Its density shape is that of a pseudo-isothermal 
ellipsoid 


Ph 


= Po[l 


(R^+Z^lq^) 


]■' 


whose characteristic parameters are also listed in Table The 
values of the central density po and the halo flattening q were 
those derived by the dynamical model from [Combes & Arn- 
[aboldi[ ( [T996l ). 

The sum of ah these mass components yields a 3D poten¬ 
tial in the whole space. In particular, we are interested in the 
resulting closed orbits in the two planes, equatorial and po¬ 
lar, which have lost their axisymmetry. In each plane, the po¬ 
tential is computed, and decomposed in Eourier components 
0(/?, 0) = Oo(/?) + O^(/?)c6>5'(m0 - 0^). We consider only the 
m = 2,4, and 6 components. 

Eohowing the epicyclic approximation, the closed orbits in 
a non-axisymmetric potential, such as a weak bar, can be devel- 
ope d in polar coordinates (/?, 0) as in [Binney & Tremam^ ( [1987| ) 
and [Schoenmakers et al.[ ( p^97] ): 


R{t) = Rq + Ri it) , 0(0 = 0o(O + 01 (0 

We take 0 = 0 as the long axis of the potential. Eor the unper¬ 
turbed circular orbit. 


0o(O = ^ 0 ^ 


/?i(0o) = Ccos(m0o), 

where m is the Eourier component in which the potential has 
been decomposed. 


C=ziA+2^) 

dR R ^^0’ 


and 


A = kI- 

with kq the unperturbed epicyclic frequency. 
The expression of the perturbed angle is 


sin(m0o) r 
01 = - 


2 


L ARq dR 




Computing derivatives of these quantities, the velocities in the 
plane of the closed orbits can be obtained as 


T/ T/ ^ o ^ 


V^ = Vo+ V„ cos(wi^i))[^( 




ARq dR 




To sum the potential of the various components, we adopt a grid 
of 480x480 in each plane, equatorial and polar. The system of 
closed orbits found in the equatorial and polar planes are dis¬ 
played in Eig.[^ 
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Fig. 7. Maps of the line-of-sight velocity (left panel) and velocity dispersion (right panel) of the gas in NGC 4650A. The cross marks the centre of 
the galaxy where rotation is V = 2880 + 19 km/s and velocity dispersion is cr = 100 + 20 km/s. In the left panel, the dashed and dotted arcs trace 
the structure of the spiral arms, as in Fig.[^ 


We then undertake to project the two components, host 
galaxy and polar disk, taking their thickness into account, with 
the adopted shapes and parameters indicated in Table As for 
the velocities, we take velocity dispersions of the stars into ac¬ 
count, as well as their asymmetric drift, i.e. we adopt for an 
axisymmetric exponential disk (with scalelength r^) the Jeans 
equation, 


= ^r(l - r/rd) - 


dcrl 

dlnr' 


where the dispersion and cr(p in the radial and tangential di¬ 
rections are related in the epicyclic approximation by 


O' 0 
CTr 


K 

2Q' 


Since the MUSE data have further confirmed that the veloc¬ 
ity dispersion profile along the major axis of the host galaxy is 
fiat with radius (see Fig. [^, we have adopted here the constant 
value crr= 65 km/s. For the z dispersion, we simply adopt 

cr^ = InGhLir), 

where S is the disk surface density. 


Since neither the host disk nor the polar disk are exactly 
edge-on, we projected them between 5° and 15° from edge-on, 
and the polar disk is not quite perpendicular to the host but in¬ 
clined by 20° from the pole. 

The model was convolved with the seeing and the final result 
is shown in Fig. and in the bottom left panel in Fig.[^ 

The main result of this analysis is that the above model can 
account for the overall 2D field of the velocity and velocity dis¬ 
persion derived by the MUSE data (see Fig. and Fig. |7]). In 
particular, the line-of-sight velocities of stars and ionized gas ob¬ 
tained by the model show the same pattern and values as those 
observed in both HG and polar disk. In the region of the cen¬ 
tral host disk, the velocity dispersion of the stars (see bottom left 
panel in Fig. [T^ is nearly constant at all radii, with comparable 
values (40 < cr < 70 km/s) to those observed in the MUSE data 
(see right panel in Fig. 

We noticed some differences between the data and the model. 
In the region of the polar disk, the velocity dispersion of the stars 
and the gas derived from the model does not completely match 
that observed in the data. The model shows a lower stellar ve¬ 
locity dispersion (~ 30 km/s, see bottom left panel in Fig. 
than that observed in the MUSE data (~ 100 - 150 km/s, see 
right panel in Fig.[^. We noticed an increase in the velocity dis¬ 
persion of the model only in a restricted region at about 3 kpc 
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Fig. 8. Left panel - Line-of-sight velocity (bottom panel) and velocity dispersion (upper panel) of the gas along the polar disk major axis, PA. = 
160°. Right panel - Zoom of the gas rotation (bottom) and dispersion (upper) curves towards the centre of the galaxy. 


(14 arcsec) on the NE and SW sides along the polar disk, where 
cr ~ 90 km/s, followed by a steep decrease. 

The velocity dispersion of the ionized gas derived from the 
model shows a gradient (see bottom right panel in Fig. [^. It 
is higher in the centre (~ 50 km/s for ^ < 2 kpc ~ 10 arcsec) 
and then it decreases at larger radii, where cr ~ 20 - 40 km/s for 
2 < R < 10 kpc (~ 10-50 arcsec). This gradient is not observed 
in the MUSE data, where the velocity dispersion of the gas along 
the polar disk is almost constant at all radii and cr ~ 40 km/s (see 
Fig. [7] right panel). 

Such differences in the velocity dispersion of stars and gas 
could arise because the real structure of the polar disk is more 
complicated than that assumed in the simple geometrical model 
described above. In particular, one of the dominant features is the 
presence of the spiral arms, which are not included in the model. 
Moreover, the data show that the polar disk is characterized by 
several star forming regions that dominate the emissions along 
the spiral arms (see Fig.[T]). Hence, the outflows from these re¬ 
gions may influence the measured velocity dispersion of the gas. 
This can explain why in the MUSE data cr has a constant value 
at ~ 40 km/s at all radii, while in the model it decreases with 
radius to ~ 20 km/s. 

Finally, in the nuclear regions (R < 2 arcsec ~ 0.4 kpc), the 
model does not show the apparent counter-rotation or the high 
velocity dispersion observed in the gas rotation curve and veloc¬ 
ity dispersion profiles (see right panel in Fig. [^. Thus, as sug¬ 
gested in Sect. |3.5[ this could be due to the presence of a real 
decoupled component or to the chaotic motions of the gas in the 
nuclear regions, also generated by an outflow, and the two fea¬ 
tures are not included in the model. Alternatively, the observed 
feature in the nuclear regions could be due to non-axisymmetric 
distortions generated by the two disks in the two planes. In fact, 
the model shows that none of the two disks have circular mo¬ 
tions, the potential is never axisymmetric in the host disk nor in 
the polar disk. This non-axisimmetry, combined with a position 
angle that is not completely perpendicular, leads to non-circular 
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motions that are seen tilted in projection (like an undetected bar 
misaligned from the major/minor axis). 

In the central regions of the HG, the model does not show 
the increasing velocity along the HG minor axis (see bottom left 
panel in Fig.[^ as observed in the MUSE data (see top left panel 
in Fig. 1^. Along the NW-SE directions, the systemic velocity 
contour in the stellar velocity field of the model is parallel to 
the HG minor axis for /? < 2 kpc (~ 10 arcsec, see middle left 
panel in Fig.[T^. Since the model does not include any rotation 
along this direction, this could be a real feature in the structure 
of the HG disk, and it does not result from the geometry of the 
two disks and its projection on the sky. This confirms the previ¬ 
ous results by |Coccato et al.| ( [2Q14| ). A definitive conclusion on 
this point will come from the 2D spectral decomposition of the 
MUSE stellar velocity field, which is the subject of a forthcom¬ 
ing paper. 


5. Summary and conclusions 

In this paper we presented new MUSE data obtained for the po¬ 
lar disk galaxy NGC 4650A. The 2D velocity distribution of 
stars and ionized gas is mapped all the way to the galaxy centre 
where the two components co-exist, overcoming the limited spa¬ 
tial coverage of the long-slit data available for this object. They 
confirm that NGC 4650A is made of two perpendicular disks that 
rotate on orthogonal planes, and that drive the kinematics right 
into the very centre of this galaxy. In the following we summa¬ 
rize the main results obtained by the 2D map of the line-of-sight 
velocity and velocity dispersion of stars and ionized gas, and we 
address the main implications on the structure and formation of 
NGC 4650A. 

- The final MUSE mosaic of NGC 4650A covers an area of 
1.5 X 2.5 arcmin (see Fig.[^. The MUSE image of the emis¬ 
sion lines (see right panel in Fig.[^ confirms that all the gas 
is associated with the polar disk. Strong emissions are ob- 
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Fig. 9. Closed orbit positions (left panels) and velocities (right) for the NGC 4650A host galaxy (top) and polar structure (bottom). 


served both in the outer arms and close to the centre of the 
galaxy. 

- Overall, the velocity field of stars inside the HG has the typ¬ 
ical pattern of a rotating disk, with receding velocities on 
the SW side and approaching velocities on the NE side. In 
addition, we also observed that i) the velocity field is not 
symmetric with respect to the centre, since the rotation curve 
is flatter on the SW side, while it is still rising on the NE 
side (see Eig.[^; and ii) there is a velocity gradient along the 
photometric minor axis of the HG. 

- The velocity dispersion remains almost constant, at (T ~ 50- 
60 km/s, at all distances from the galaxy centre and RA.s 
inside the HG, while it increases along the polar disk, up to 
(T ~ 100 km/s in the outer regions. 

- Inside the HG regions, the photometric and kinematic RA. 
differ by about 5 degrees (see Eig. [^, with the photometric 
R.A., PA.pho = 67 + 2 degrees, is smaller than the kinematic 
RA. 

- The 2D line-of-sight velocity of the ionized gas, which is 
distributed along the polar disk of NGC 4650A, (Eig.[^ left 
panel) is mapped up to 75 arcsec (~ 16 kpc) from the cen¬ 


tre. It shows a large amount of differential rotation (of about 
100- 120 km/s) from the centre up to the outer regions. On 
the NW and SE sides of the polar disk, the velocity of the 
approaching and receding arms is measured, including the 
regions where they overlap with the underlying disk. 

- The gas velocity dispersion remains almost constant (cTgas ~ 
30 - 40 km/s) at all radii (see right panel in Eig. [7]), except 
inside 2 arcsec from the galaxy centre, where we measured a 
steep increment of cr that reaches a value of about 110 km/s. 
In the same range of radii, we also observe a hint of a de¬ 
coupled component in the velocity curve (see right panel in 
Fig.§. 


In order to reproduce the observed structure and the new 2D 
kinematics for NGC 4650A given by the MUSE data, we con¬ 
structed a multicomponent mass model made by the combined 
projection of two disks, plus a dark matter hal o (see Sect. [4]). It 


is based on th e dynamical model developed by [Combes & Arn 
|aboldi| ( p~996| ), where the dark halo is flattened and its major axis 
is aligned with the polar disk. The structural parameters for the 
light distribution take the extended and deeper photometry pub- 
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Fig. 10. Projection of the model of the stellar (left panels) and ionized gas kinematics (right) for NGC 4650A. The surface density distribution 
of stars and gas are shown in the upper panels, the projected velocity in the middle panels, and the velocity dispersion in the bottom panels. The 
colour scales are indicated in km/s. North is up, east is on the left. 
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Table 1. Masses and scales for the model components 


Component 

M 

h 

r 

n 

ri 


10® Mo 

[kpc] 

[kpc] 

[kpc] 

[kpc] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

HG bulge 

Mb 


rB 




0.2 


0.17 



HG disk 

Md 

hd 

rd 




10.3 

0.5 

0.948 



Polar disk 

MpD 

hpD 





15. 

0.5 


5.95 

6.8 

HI disk 

Mhi 

hm 





7.2 

0.5 


3.4 

15.3 

DM halo 

Mdh 

z 

rh 




15 

1.2 

6.0 




Notes. Col. 1: Component included in the mass model. Col. 2: To¬ 
tal mass of each component. Col. 3: Characteristic scale-height of 
each component. Col. 4: Characteristic scale-length of each component. 
Col. 5 and Col. 6: ri and r 2 are the internal and external radii of the stel¬ 
lar and gaseous polar disks. 


lished by |Iodice et aL] ( |2002b| ) into account. It is a “minimum” 
model that does not include the spiral arms in the polar disk, the 
rotation along the minor axis of the HG and the nuclear decou¬ 
pled component in the centre. 

By comparing the observations with the 2D kinematics de¬ 
rived by the model, we find that the complex velocity field re¬ 
vealed by the MUSE data for NGC 4650A is overall well repro¬ 
duced by this simple model, also in the central regions of the 
galaxy, where the two components coexist (see Fig.p^. 

This result is a strong constraint on the dynamics and forma¬ 
tion history of the polar disk, in particular on the nature of the 
dark matter. In the framework of disk galaxy formation, an ex¬ 
tended and massive polar disk like that observed in NGC 4650A 
can form through the accretion of cold gas along a filament. In 
this scenario, the central HG is disk that formed first. The decou¬ 
pling of the angular momentum at a certain epoch of the forma¬ 
tion history moved the accretion of gas in the orthogonal direc¬ 
tion, forming the polar disk ( [Snaith et al.|20I2t|Combes|20I4| ). 
Alternatively, a polar disk could also form from the tidal accre¬ 
tion of material from a gas-rich companion galaxy. Both mech¬ 
anisms were successfully tested for NGC 4650A (|Spavone et al 

[MTOl ). 


The analysis performed in this work shows that a mass model 
made by the two perpendicular disks reproduce well the ob¬ 
served structure and kinematics of NGC 4650A, which result 
from one of the accretion mechanism described above. Since the 
mass model includes a fiattened dark halo, having its major axis 
aligned with the polar disk, this suggests that if the polar disk is 
formed through the accretion of material from outside, dissipa¬ 
tive dark matter might to be essential. 


Some differences are noted between the observed and mod¬ 
elled kinematics. They are the rotation velocity along the minor 
axis of the HG, the high stellar velocity dispersion in the polar 
disk, and a possible decoupled component in the galaxy centre, 
which are detected in the data but not in the model. 


They do not affect the main conclusion discussed above, 
since they are limited to regions that are not dynamically im¬ 
portant for the whole model. In fact, inside the HG, the minor 
axis rotation is for /? < 5 arcsec (~ 1 kpc), which are only about 


15% of the measured kinematics along the major axis, and which 
are fully reproduced by the model. 

In the polar disk, we measured an increasing velocity disper¬ 
sion of stars, for 5 < /? < 10 arcsec (~ 1 - 2 kpc), out to about 
110 km/s (see right panel in Fig.|^, and a difference in the ve¬ 
locity dispersion of the gas of about 20 km/s at 16 kpc (see right 
panel in Fig.[7]and right bottom panel in Fig.|^. Along the polar 
disk, the dynamics is strongly constrained by the high velocities 
(~ 120 km/s) of the HI gas at 40 kpc from the centre. 

As pointed out in Sect.|^ such differences in both stellar and 
gaseous velocity dispersion could be a consequence of the per¬ 
turbed structure of the polar disk, which are not taken into ac¬ 
count by the simple model made in this work. The differences in 
the measured velocity dispersion of the gas with respect to that 
derived from the model could be due to the outflows from the 
several star forming regions that dominate the emissions along 
the spiral arms. On the other hand, the higher velocity disper¬ 
sion of the stars in the outer regions of the polar disk and the 
counter-rotation in the centre of the polar disk, if the latter is a 
real decoupled component, could be the result of the angular mo¬ 
mentum evolution of the accreted material during the formation 
process. 
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